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Abstract 

Recent strong lensing data and rotation curves of dwarf galaxies indicate that 
many galactic clusters may have a soft core instead of a central cusp in their 
density distribution. This result challenges the standard CDM (Cold Dark 
Matter) theory prediction based on N-body simulations. We find that the 
observed density profile is consistent with that of a spherical gas of degenerate 
fermions such as neutrinos in hydrostatic equilibrium. Also, we compare two 
different models of the dark matter halo and their predictions for the hot gas 
profile. 
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Recent strong lensing data and rotation curves of dwarf galaxies suggest that matter 
distribution takes the form of a soft core at the centers of some galaxy clusters [1,2]. How- 
ever, N-body simulations based on the CDM theory predict that the density profile of the 
dark matter halo should be singular at the center which is called the NFW profile [3,4]. 
One possibility to resolve the discrepancy is that the dark matter may not be coUisionless 
as assumed in previous simulations. They may interact with each other and form a profile 
different from the NFW profile. In the following, we show that indeed the observed den- 
sity profile is consistent with that of a degenerate gas of massive neutrinos in hydrostatic 
equilibrium. We will first review a method of inferring the dark matter distribution from 
the density profile of the hot gas in a galaxy cluster [5,6]. Then we will demonstrate that 
one of the dark matter candidates, massive neutrinos, which can be either coUisionless or 
self-interacting, can account for the observed profiles. 

One of the most common methods to probe the dark matter distribution is to make use 

of the gas density profile. We assume that the hot gas is in hydrostatic equilibrium, i.e., 

dng{r) _ GM{r)mgng{r) 
dr ~ kTr^ ' ^ ^ 

where ng{r) is the number density of hot gas, M(r) is the total mass enclosed within a radius 

r, k is the Boltzmann constant, and nig is the mass of the gas particles. Here we assume that 

the hot gas is isothermal throughout the cluster, with a constant temperature T, which is 

valid as the heat conduction by electrons is quite efficient. By knowing the number density 

profile of the hot gas, we can probe the distribution of dark matter. There is a general form 

of the hot gas density profile called the P model [5], 

where f = r/rc, Tc being the core radius, /3 is a parameter, and Uq is the central number 
density. Different clusters have different values of no, Vc and /3. Now, we can solve Eq. (1) 
and Eq. (2) and get: 



For r much greater than Tc, M(r) oc r, whereas for r much smaller than Tc, M(r) oc r^, 
which means that the density is constant, 

In general, from Eq. (3), we can obtain the mass density of the dark matter, assuming that 
the total mass is dominated by dark matter together with the hot gas, 

which is quite different from the NFW profile, 

Pdm = / / w|", — , (6) 

where po and are parameters. Possible sources of the discrepancy are problems in the 
NFW profile, or the /3 model, or the assumption of isothermal distribution in the cluster. 
For example, it has been proposed that including the cooling flow of hot gas in the above 
model will result in a dark matter distribution that is closed to the NFW prediction [7]. 

Nevertheless, the hot-gas method is model dependent and the uncertainties are still large. 
To understand more about the dark matter, one powerful and model independent alternative 
is to use strong lensing to probe the distribution of dark matter. This method still has large 
errors and is limited in data collection. One of the clusters studied recently is CL0024-I-1654 
[1]. It was concluded that CL0024-I-1654 has a soft core in the center which is different from 
the NFW profile [8]. Furthermore, some observations of dwarf galaxies mass distribution 
also indicate that the inner halo density is nearly constant [2]. 

It has been proposed that cold dark matter may have weak self-interaction instead of be- 
ing coUisionless [9]. Following this suggestion, we assume that the dark matter can establish 
hydrostatic equilibrium, 

ciPpM ^ G'M(r)pDM(r) 
dr ' 

where Pqm is the dark matter pressure. Then, we can combine Eq. (1) and Eq. (7) to get 



(7) 



^-Pdm(pdm) _ kT^dUg 

PDM Ug 

Prom Eq. (8), we can obtain the equation of state (EOS) of dark matter Pdm(pdm) from the 
hot gas density profile. If we assume a polytropic relation between the pressure and density 
of the dark matter, 

PdM = i^PDM^ , (9) 

where K and 7 are parameters in the EOS, Eq. (8) becomes 

(10) 



Pdm — PDMc 



,PgcJ kT{-f - 1) 

with 7 > 1, and pgc and pdmc are central density of hot gas and dark matter respectively. 
The hot gas density profile directly tracks the dark matter profile and is sensitive to the 
steepness of the dark matter EOS. 

In the following, we consider massive neutrinos as a specific model of dark matter and 
solve the cluster core problem. At low density, neutrinos are coUisionless as they have very 
small interaction cross-section. However, at high density they become degenerate, and they 
can therefore be in hydrostatic equilibrium. It has been proposed that such neutrino "stars" 
may be found in the centers of galaxies [10,11]. Here we propose that they may exist at the 
centers of galaxy clusters. Since the interaction cross-section of neutrinos is small, we can 
treat them as a zero temperature gas. The degenerate pressure for non-relativistic neutrinos 
is given by 

' , (11) 



where rii, and rrii, are the number density and mass of the neutrinos. We can then solve for 
the density profile of a neutrino star assuming hydrostatic equilibrium (shown in Fig. 1), 
which is expectedly similar to that of a white dwarf. Here we set the neutrino mass to be 
3 eV and the central number density to be 10^ cm~^, which fits the observational data for 
clusters. The radius of the neutrino star is several kpc which is comparable to the cluster 



scale. The total mass is about 10^^ solar masses, which is also close to the total mass of a 
cluster. A neutrino star has a core in the center, the radius of which depends on the central 
mass density only. The higher the central mass density is, the smaller is the core radius 
and vice versa. The inner density profile is similar to the data obtained from gravitational 
lensing. This suggests that at least some clusters may be degenerate fermions dominated. 
From Fig. 1, we can approximate the neutrino star density profile as following: for r << ro, 
M(r) = Cr^, where C is a constant, and for r >> ro, M(r) = Mq, where Mq and ro are 
the total mass and core radius of the neutrino star. Using the above approximation, we can 
analyse the density profile of the hot gas inside a cluster core if its mass is dominated by 
the neutrino star. If M(r) = Cr^, we have 

n^^noe-"'^', (12) 

where a — GCmg/2kT. For r << a~^/^, the hot gas profile is smooth and flat. Therefore, 
the core radius is 

Clearly, rc is not necessarily equal to ro, but it is affected by the central mass density of the 
neutrino star and temperature of the hot gas. In our model, for T = 10® K and central mass 
density Pc — 10~^^ gcm~^, Vc is about 150 kpc. 

We can also calculate the density profile of the hot gas numerically. In Fig. 2, we show 
the hot gas profile using the above parameters of a neutrino star, and we find that it is 
well fitted by using a /3 model with /? = 1.1 and rc=150 kpc. This indeed agrees with the 
observational value of P between 0.9-1.1 [12,13]. 

We now compare the results of using three different models, the NFW dark matter 
profile, neutrino star model and hot gas self-bounded profile. Here, we continue to assume 
that the temperature of the hot gas is constant throughout the cluster. Fig. 3 shows the 
hot gas profiles from these three models. We can see that the three hot gas profiles have 
similar central mass density and temperature. However, the slopes of the hot gas profiles 



for r > Tc are different, being smallest for self-bounded model and largest for the neutrino 
star model. In fact, the slope of the hot gas profile depends not only on the dark matter 
profile but also on the temperature profile of the hot gas and velocity dispersion of the dark 
matter particles. Therefore, the dark matter profile obtained by using the hot gas profile 
alone is model dependent. For r > tq, the total density profile is dominated by the hot gas 
and the total density tail pt is approximately proportional to which agrees with the 
observational result [14]. To understand the dark matter distribution better, gravitational 
lensing is a potentially more powerful method. It is model independent, and we can know 
the mass density profile directly without using any hot gas model if the cluster is dominated 
by the dark matter. For example, the cluster CL0024-I-1654 observed by Tyson et al. does 
not have a cD galaxy, which means that the mass of dark matter is dominant at the center 
of the cluster [2]. Gravitational lensing data suggests that the NFW profile may not be 
correct. 

We have shown in this article that observational data on cluster hot gas is consistent with 
the existence of neutrino stars, which can account for at least parts of the dark matter. But 
how are these neutrino stars formed? It is commonly believed that there is a cosmological 
neutrino background originating from the Big Bang [15], corresponding to a temperature 
of about 1 K now, or about 0.1 meV in energy scale. Therefore it is conceivable that 
cosmological neutrinos with rest mass even in the sub-eV range are non-relativistic. Recent 
experimental observations of neutrino oscillations point to a range of the mass difference 
between muon and tau neutrinos to be 5 x 10^"^ < Am^ < 6 x 10^^ eV^ [16]. Therefore at 
least one species of the cosmologocal background neutrinos should be non-relativistic, and 
the interplay between gravitational attraction and degenerate pressure makes it possible to 
create neutrino stars. We have carried out hydrodynamics simulations of their formation 
process, and preliminary results indicate that hydrostatic neutrino stars can indeed form 
with a wide range of parameters [17]. 

One of the interesting properties of the neutrino star is that it has a fiat core, which 
agrees with the gravitational lensing data as well as the rotation curves of dwarf galaxies, 
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but contradicts with the NFW profile. Recently, Sand et al. has provided an upper limit 
of the slope of the central density profile a < 0.57 which is different from NFW's (a = 1) 
[18]. In our model, a is approximately zero. In order to solve the "core problem" of galactic 
clusters, more data will clearly be needed. 

The work described in this paper was substantially supported by a grant from the Re- 
search Grants Council of the Hong Kong Special Administrative Region, China (Project No. 
400803). 
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FIG. 1. The mass density profile of a neutrino star with central mass density = 10^^^ gcm~^ 
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FIG. 2. Numerical result of the hot gas profile and beta model fitting with T=10^ K, rc=150 kpc 
and 71^=20000 at the center. Right panel: hot gas profile together with beta model fitting. Left 
panel: the central core of the hot eras in loe scale: darl<: matter profile is also shown (solid line). 
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FIG. 3. Hot gas profiles in four different models: NFW model from CDM simulation, with 
rs=300 kpc (dot-dashed line), neutrino star model (solid line), beta model (dashed line), and 
self-bounded model (dotted line). For the latter three models, same parameters as in Fig. 2 have 
beer "^^^^ 
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FIG. 4. Dark matter profiles in neutrino star model (solid line) vs. NFW model from CDM 
simulation, with rs=300 kpc (dot-dashed line). Same parameters as in Fig. 2 have been used. 
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